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ABSTRACT 

Based on the standard energy reservoir of gamma-ray burst (GRB) jets and the re- 
lationship between peak energy (Ep) of vF^ spectrum and the equivalent-isotropic 
energy (Eiso), we test the GRB probabiUty distributions as a function of viewing an- 
gle {9) and of viewing angle together with redshift z, P{9) and P{9, z), predicted by 
the universal-structured jet (USJ) model with simulations on observational bases for 
a detection threshold of S" > 4 x 10~^ erg cm^^, where S is the total fluence in an 
energy band > 20 keV. We simulate a sample including 10^ GRBs, and then derive 
P{9) and P{6,z). We show that both simulated and theoretical P{9) and P{9,z) are 
consistent with each other. This consistency seems to be regarded as support for the 
USJ model. A small difference between the theoretical and simulated results might be 
due to the observed Ep distribution used in our simulations, which perhaps suffers a 
httle of biases for E'p < 100 keV. 

Key words: gamma rays: bursts — gamma rays: observations — ISM: jets and 
outflows — methods: statistical 



1 INTRODUCTION 

Significant progress on understanding the nature of gamma- 
ray bursts (GRBs) and their afterglows has been made in 
the recent decade (for reviews see Fishman & Meegan 1995; 
Piran 1999; van Paradijs et al. 2000; Cheng & Lu 2001; 
Meszaros 2002). However, this phenomenon is still of a great 
mystery. One of the most difficult parts of solving this mys- 
tery is their central engines, which are far away from emis- 
sion regions and are hid behind our observations (e.g., Cheng 
& Lu 2001). It is believed that GRBs are produced by conical 
ejecta powered by the central engines. Sharp breaks and/or 
quick decays in afterglow light curves are thought to be ev- 
idence for coUimated ejecta (e.g., Rhoads 1999; Sari, Piran, 
& Halpern 1999; Dai & Cheng 2001). Starting with GRB 
990123 such evidence has been rapidly growing up (e.g., 
Bloom, Frail, & Kulkarni 2003 and references therein). In 
a previous work we showed that the current observations of 
prompt gamma-ray emissions and X-ray afterglow emissions 
satisfy the prediction of GRB jet (Liang 2004). Both the 
dynamics and afterglow light curve behaviors of GRB jets 
were investigated numerically (Panaitescu, Meszaros & Rees 
1998; Huang et al. 2000; Moderski, Sikora & Bulik 2000; 
Granot & Kumar 2003; Kumar and Granot 2003; Salmonson 
2003), and some GRB afterglows of interest were also fitted 
(Panaitescu & Kumar 2001a, 2001b, 2002). 

GRB jet structure models are currently under heavy 



debate. Two currently competing models are universal- 
structured jet (USJ) model (e.g., Meszaros, Rees & Wijers 
1998; Dai & Gou 2001; Rossi, Lazzati, & Rees 2002; Zhang 
& Meszaros 2002a; Granot & Kumar 2003; Kumar & Gra- 
not 2003; Panaitescu & Kumar 2003; Wei & Jin 2003) and 
uniform jet model (e.g., Rhoads 1999; Frail et al. 2001). 
X-ray flashes (XRFs), which are thought to be the low- 
energy extension of typical GRBs (Kippen et al. 2003; Lamb, 
Donaghy, & Graziani 2003; Sakamoto et al. 2004), seem to 
provide more jet structure signatures. If both GRBs and 
XRFs are the same phenomenon, any jet structure model 
should present a unified description for GRBs and XRFs. 
Zhang et al. (2004a) showed that the current GRB/XRF 
prompt emission/afterglow data can be described by a quasi- 
universal Gaussian-like (or similar structure) structured jet 
with a typical opening angle of ~ 6° and with a standard 
jet energy of ~ 10^^ ergs. Based on High Energy Tran- 
sient Explorer 2 {HETE-2) observations. Lamb et al. (2003) 
proposed that the uniform jet model reasonably describes 
the unified scheme of GRBs/XRFs. Very recently, the two- 
component jet model was advocated by Berger et al. (2003a) 
based on the observations of GRB 030329, which has two 
different jet breaks in an early optical afterglow light curve 
(0.55 days; Price et al. 2003) and in a late radio light curve 
(9.8 days; Berger et al. 2003a). Millimeter observations fur- 
ther support the two-component jet in this burst (Sheth 
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et al. 2003). Simulations of the propagation and eruption 
of relativistic jets in massive Wolf-Rayet stars by Zhang, 
Woosley, & Hoger (2004b) also showed signatures of the 
two-component jet. In a previous work, we found a bimodal 
distribution of the observed peak energy (Ep) of i^F,^ spec- 
tra of GRBs/XRFs, and concluded that the two-component 
jet model can explain this distribution (Liang &i Dai 2004). 
Our results suggest that the two-component jet seems to be 
universal for GRBs/XRFs. Numerical calculations of such a 
model were also presented by Huang et al. (2004). 

The USJ model has a predictive power. This might be 
served as a test for this model. Perna, Sari, & Prail (2003) 
found that both the observational and theoretical probabil- 
ity distributions of 6, P(0), are roughly consistent with each 
other. However, Nakar, Granot, & Guetta (2004) showed 
that the theoretically 2-dimensional probability distribution, 
P{z,6), poorly agrees with the observational one. We note 
that both Perna et al. (2003) and Nakar et al. (2004) came 
to conclusions from the comparison of the USJ model's pre- 
dictions with the current observations. Whereas, the current 
sample of GRBs with available is very small (16 GRBs in 
their works). Observational biases and sample incomplete- 
ness affect greatly this sample. It is difficult to draw a ro- 
bust conclusion in a statistical sense from the comparison 
of P{6) or P{9, z) derived from this sample with that pre- 
dicted by the USJ model. In addition, both theoretical P{d) 
and P{9, z) are sensitive to instrument's threshold setting. 
In Perna et al. (2003) and Nakar et al. (2004), the thresh- 
old is the same as CRGO/BATSE. However, most of the 
bursts with 6 available in the current GRB sample were 
not observed by CGRO/BATSE. This difference also affects 
greatly the comparison results. We therefore make a further 
test by simulations. Based on the standard energy reservoir 
of GRB jets (Frail et al. 2001; Bloom et al. 2003; Panaitescu 
& Kumar 2001a; Piran et al. 2001; Berger, Kulkarni, & 
Frail 2003b) and the relationship between Ep and -Eiso, 
the equivalent-isotropic energy (Amati et al. 2002; Lloyd- 
Ronning & Ramirez-Ruiz 2002; Atteia 2003; Sakamoto et 
al. 2004; Lamb et al. 2003; Liang, Dai , & Wu 2004; Yone- 
toku et al. 2004), we test the P{e) and P(6», z) predicted by 
the USJ model with simulations on observational bases for 
a detection threshold of S* > 4 x 10~^ erg cm~^, where S is 
the total fluence in an energy band > 20 keV. We show that 
both P{9) and P{9, z) derived from simulations and from 
the USJ model are consistent with each other. 

This paper is arranged as follows. The empirical model 
is presented in section 2. Our simulation procedure is de- 
scribed in section 3. The results of the comparison of our 
simulation results with that predicted by the USJ model are 
presented in sections 4. Discussion and conclusion are pre- 
sented in section 5. Throughout this work we adopt Hq = 65 
km s"^ Mpc -\ = 0.3, and Oa = 0.7. 



2 MODEL 

Our simulation analysis model is based on the statistical 
results of a standard energy reservoir in GRB jets and the 
relationship between Ep and Eieo- The energy release of a 
GRB jet in the gamma-ray band is given by 

£52 = £^iso,52(l - C0S9), (1) 



where i?iso,52 = -Biso/lO^^ergs. The relationship between 
Eiso and E-p is given by 

Eiso,52 = k[Ep,2{l + z)f , (2) 

where fc is a coefficient, and E-p^-z — -Ep/lO'^keV. 
This empirical relation was further investigated in 
standard synchrotron/inverse-Compton/synchro-Compton 
models (Zhang & Meszaxos 2002b; Ramirez-Ruiz & Lloyd- 
Ronning 2002). Sakamoto et al. (2004) and Lamb et al. 
(2003) pointed out that HETE-2 observations not only con- 
firm this correlation, but also extend it to XRFs. In a recent 
work, we also proved that this relationship holds within a 
burst (Liang, Dai, & Wu 2004). From Eqs. (1) and (2), we 
obtain 

. = arccos{l- ^^^J-^^^jJ . (3) 

Since the values of k and £52 can be obtained by the current 
GRB sample, the 9 value of a burst can be estimated by Eq. 
(3), if its Ep and z are available. If both distributions of Ep 
and z can be established by observations or theory, one can 
perform simulations to derive a GRB sample at the least 
observational bias, and obtain P{9) and P{9, z) from the 
simulated sample. Such a Ep distribution for a given detec- 
tion threshold may be established by the BATSE observer 
tions. The redshift distribution is impossible to be modelled 
by the presented GRB sample because of various observa- 
tional biases, such as GRB detection, localization, and red- 
shift measurement (Bloom 2003). Please note that, in our 
analysis, we focus on the comparison of the simulated P{()) 
and P(S, z) with that predicted by the USJ model. If we 
use the same redshift distribution model in our simulations 
and in the USJ model, the comparison results are not signif- 
icantly affected by the redshift distribution (for a detailed 
discussion in section 5). Thus, we perform simulations to de- 
rive P{6) and P{9, z) and compare them to the predictions 
of the USJ model for a given detection threshold, assuming 
the GRB rate as a function of redshift is proportional to the 
star formation rate. We present the distributions of these 
parameters used in our simulations in the following. 

The strong correlations shown in Eqs. (1) and (2) in- 
dicate that both £52 and k have a standard candle feature. 
They should not be significantly affected by observational bi- 
ases and the sample selection effect. Their distributions can 
be modelled by the present GRB sample. ^From a sample 
of GRBs with 9 available presented by Bloom et al. (2003), 
we find that the £52 distribution can be modelled by a log- 
normal one with a center at —0.9 and a criogg^j = 0.3, i.e., 

0°g£52+2_9)i 

P(Ioge52) = Pe,oe 05 (4) 

where Pe,o is a normalized constant given by 
P(log£52)dlog£52 = 1. Prom a GRB sample presented 
by Amati et al. (2002), we find that the k distribution can 
be modelled by 

P(fc) = Pk,oe-'^^, (5) 

where Pk,o is a normalized constant given by Jp°° P{k)dk = 
1. 

The observed Ep distribution is significantly affected by 
observational biases and sample selection effects, especially 
when the completeness at low fiuxes is considered. We mod- 
elled the Ep distribution for a given detection threshold from 
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GRB observations. In a previous paper, we studied the ob- 
served Ep distribution of GRBs and XRFs, combined with 
both HETE-2 and BATSE observations. We found that the 
observed Ep distribution for GRBs/XRFs is a bimodal one 
with peaks of ~ 30 kcV and - 200 kcV (Liang & Dai 2004). 
The ~ 30 keV peak has a sharp cutoff at its lower energy side 
because of the limitation of the current instrument thresh- 
olds. This bimodal distribution should be further examined 
by future observations. In addition, for some extremely soft 
XRFs, they violate the model used in our analysis (Eq. 3; 
e.g., Liang & Dai 2004). We therefore cannot use this Ep dis- 
tribution in our analysis. It is well known that the observed 
Ep distribution of typical GRBs is narrowly clustered at 
200-400 keV (e.g., Preece et al. 2000). The Ep distribution 
of bright GRBs has been well established. Thus, we consider 
only bright GRBs in our analysis. Preece et al. (2000) pre- 
sented the observed Ep distribution for a long, bright GRB 
sample selected by 5 > 4 x 10~^ erg cm ~^ in an energy 
band > 20 keV. We use this GRB sample to model the Ep 
distribution, which is 



P{\ogEp,2) = PE,oe' 



, (log Ep, 2-0-38)^ 



(6) 



Pe,o is a normalized constant, which is given by 
j;^^PilogEp,2)dlogEp,2 = 1. 

The GRB rate as a function of redshift is assumed to 
be proportional to the star formation rate (e.g., Bromm & 
Loeb 2002 ). The SFl model from Porciani & Madau (2001) 
is used in our analysis. Thus, the redshift distribution is 
modelled by 



+ 45 



(7) 



where Pz,o is a normalized constant. Since the largest red- 
shift of GRBs is 4.5, we restrict z < 4.5, and Pz,o is given 
by J^ 'P{z)dz = l. 



3 SIMULATION PROCEDURE 

Based on the empirical model and its parameter distribu- 
tions presented above, we simulate a sample including 10^ 
GRBs. Each simulated burst is characterized by a set of 
(fc, £52, Ep^2, z), and it also satisfies our threshold setting. 
The threshold setting corresponds to the sample selection 
criterion of the sample in Preece et al. (2000), which is 
S > 4 X 10~^ erg cm~^. Our simulation procedure is de- 
scribed as follows. 

(1) Derive the accumulative probability distributions of 
the parameters from Eqs. (4)-(7), Q{x), where x is one of 
these parameters. For saving calculation time, we use the 
discrete forms of these distributions. By giving a bin size 
and truncating a; at a given value a, we calculate P{xi) by 
Eqs. (4)-(7), where Xi is the ith bin of parameter x, and then 
derive Q{xi) = J]-i^q-P{x-i). The normalized constants in 
Eqs (4)- (7) arc given by Q(xi)|a;,=a = 1- 

(2) Simulate a GRB. We first generate a random num- 
ber m [0 < m < 1), and obtain the value of x from the 
inverse function of Q{x) = m, i.e., x = Q^^{rn). Please note 
that Q{x) is in a discrete form. The x value for a given 
m is restricted by Q{xi) < m and Q{xi+i) > m. Then, 
the x value is given by a; = {xi+i + Xi)/2. Repeating this 



step for each parameter, we get a simulated GRB charac- 
terizing by a set of {k, £52, Ep^2, z). If this burst violates 
e52/fe[Bp,2(l + z)]^ < 1, it does not satisfy Eq. (3) and is 
ruled out without further consideration. 

(3) Examine whether or not the simulated GRB satisfies 
our threshold setting. For a simulated GRB, we calculate its 
flux in the observer frame by F = iJiso(l + z)/A'kD\{z)T, 
where Dl{z) is the luminosity distance at z, T is the "effec- 
tive" duration of the burst, and -Eiso is given by Eq. (2). 
The "effective" duration in Perna et al. (2003) is 8 sec- 
onds. In our analysis we take T = 12 seconds (for a de- 
tail discussion in section 4). The limit of F is given by 
-Fiim = S/T = 3.3 X 10-" erg cm^^ s^^ U F > Fn^, the 
burst is included into our mock GRB sample. 

(4) Repeat steps (2) and (3) to obtain a mock sample 
of 10*^ GRBs. 



4 RESULTS 

We calculate the d value for each mock GRB by Eq. (3) . The 
relative probability of 6 is shown in Figure 1 (the step line) . 
For comparison, the results derived from the USJ model 
is also plotted in Figure 1 (the straight line). Please note 
that, when we calculate the 9 distribution by the USJ model, 
the "effective" duration T and the threshold setting axe dif- 
ferent from that used in Perna et al. (2003). T is an ad- 
justable parameter, which is taken as 8 seconds in order to 
obtain a distribution that is agreement with the observa- 
tional data in Perna et al. (2003). From Figure 1, we can see 
that the simulated distribution is a log-normal one centering 
at log 6 = —1.3. This result is significantly different from the 
observational result^-The mean of 9 for a sample of GRBs 
with 6 available in Bloom et al. (2003) is 0.12 rad, while 
our simulation result is clustered at 0.04 rad. Hence, we do 
not take the same T value as that in Perna et al. (2003). 
We adjust its value to let both the simulated 6 distribution 
and that derived by the USJ model peak at the same posi- 
tion. We obtain T ^ 12 seconds. For the threshold setting, 
we adopt Fum = 3.3 x 10"^, as we mentioned in the above 
section. The maximum redshift (zmax) up to which a burst 
with apparent angle 9 can satisfy the detection threshold is 
obtained from 



Flim = 



)(1 - cos6')T 



(8) 



where Ej = 1.33 x 10^^ ergs (Bloom et al. 2003). Since the 
maximum redshift in current GRB observations is 4.5, we 
set «max = 4.5 when «max > 4.5. 

^From Figure 1, we find that the 9 distribution of the 
mock GRB sample ranges in 0.001-0.25 rad, and the two dis- 
tributions are well consistent when 6 < 0.10 rad, while the 
simulated 9 distribution is lower than that predicted by the 
USJ model when 9 > 0.10 rad. We evaluate the consistency 
of these two distributions by a Kolmogorov-Smirnoff (K-S) 
test (Press et al. 1997, p. 617). The result of the K-S test 
is described by a statistic of Pk-s: a small value of Pk-s 

^ Please note that the simulated distribution is dependent on 
the star formation rate model and the detection threshold. It can- 
not be regarded as a true distribution. Any comparison between 
this distribution with observational results is meaningless. 
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Figure 1. The comparison of tiie relative probability distribu- 
tions derived from our simulations (the step line) with that pre- 
dicted by the US J model (the straight line). 

indicates a significant difference between two distributions 
(Pk-s = 1 indicates that two distributions are identical, and 
Pk-s < 10~* suggests that the consistency of two distribu- 
tions should be rejected; e.g., Bloom 2003). The K-S test 
shows that Pk-s = 0.73, indicating that the consistency of 
the two distribTitions is acceptable. 

Nakar ot al. (2004) found that the 2-diniensional dis- 
tribution of the probability from observations, P{\ogd, z), 
is poorly consistent with that predicted by the USJ model. 
We have mentioned in section 1 that this discrepancy might 
be due to the threshold sotting and sample incompleteness. 
We compare the 2-dimensional distribution derived from our 
simulations (the gray contour plot) with that predicted by 
the USJ model (the line contour plot) in Figure 2. The 
labels marked in Figure 2 indicate the probabilities of a 
burst within the corresponding regions predicted by the USJ 
model, and the filled colors from light gray to dark gray rep- 
resent the same probability regions derived from our simu- 
lations. From Figure 2, we find that our simulation result is 
consistent with that predicted by the USJ model. To clearly 
illustrate this consistency, we also compare the la region 
(the probability of a burst in this region is ~ 0.68) of the 
two distributions in Figure 3. It is found that the two dis- 
tributions axe most overlap. The results shown in Figures 2 
and 3 well agree with that shown in Figure 1. 



5 DISCUSSION AND CONCLUSIONS 

Based on the standard energy reservoir of GRB jets and the 

relationship between Ep and -Eiso, we test P{9) and P{6,z) 
predicted by the USJ model with observational-based simu- 
lations for a detection threshold of 5 > 4 x 10^'' erg cm~^. 
We simulate a sample including 10^ GRBs, and then derive 
P{9) and P{0, z). We find that both simulated and theoret- 
ical P{0) and P(0, z) are consistent with each other. 

As can be seen in Figures 1-3, a small difference between 
simulated and theoretical P{0) and P{0, z) results from the 
range of ^ > 0.1 rad. This difference increases as 6 increases 



Figure 2. The comparison of the 2-dimensional distribution of 
the relative probability derived from our simulations (the gray 
contour) with that predicted by the USJ model (the line contour 
)• 
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Figure 3. The comparison of the la region of the relative proba- 
bility distribution derived from our simulations (the gray contour) 
with that predicted by the USJ model (the line contour ). 

when 6 > 0.1 rad. The accumulative probability of 6 < 0.25 

rad derived from our simulations is ~ 99%, while it is 0.83 
predicted by the USJ model. From Eq. (3) one can see that 
a burst with lower Ep tends to have a larger 0. Thus, we 
suspect that the difference might be due to the limitation of 
our simulations and the observational biases of Ep distribu- 
tion. In our simulations, we derive 6 for each burst by Eq. 
(3), which requires Ep{l + z) > 100(e52/fc)"-''' ~ 32(£5i/A:)°'^ 
keV, where (e5i/A:)° '^ ~ 1. Those bursts with Ep{l + z) < 32 
keV is excluded from our simulated GRB sample. However, 
this limitation could not significantly affect our results since 
the Ep distribution used in our simulations mainly ranges 
in Ep > 100 keV. The other possible reason might be the 
observational biases of the observed Ep distribution used in 
our simulations. For a burst with 9 > 0.25 rad, its Ep value 
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is less than 100 kcV, assuming that £52 — 0.133, k = 1, 
and z = 1. Thus, the difference between the simulated and 
theoretical results at ^ > 0.25 rad might be the bias of Ep 
distribution at Ep < 100 keV. The distribution is taken 
from BATSE Largo Area Dotoctor (LAD) observations. The 
BATSE/LAD thresholds are sensitive to the spectral shape 
of a burst: the threshold for a hard burst is lower than the 
threshold for a softer burst Bursts with low Ep seems to be 
more easily lost by BATSE/LADs. The BATSE/LADs op- 
erate in the 20-2000 keV band but usually trigger in the 50- 
300 keV band (e.g., Band 2003). BATSE generates a GRB 
record when the observed photon fluxes of 2 or 3 energy 
bands (normally 50-100 kcV and 100-300 KeV bands) arc 
simultaneously beyond the thresholds. In our simulations, 
we take a threshold as Fth = S/T = 3.3 x 10"® erg cm"^ 
s~^ . We examine whether or not this flux is beyond the LAD 
thresholds. We assume that the GRB spectral shape is the 
Band function (Band 2003) with a low-energy index a = — 1 
and a high-energy index f3 — —2, and then translate the Fth 
in units of erg cm~'^ s~^ to photon cm~^ s~^ in a timescale 
of 1.024 seconds. We obtain 4.57 and 3.29 photon cm~^ s~^ 
in the 50-100 keV and 100-300 keV bands when Ep = 100 
keV, respectively. A normal GRB with such a photon flux is 
intense enough triggering the BATSE/LADs. However, The 
observed Ep distribution by BATSE/LADs mainly ranges 
in 100 - 1000 keV. This might suggest that BATSE/LAD 
sensitivity decreases significantly for Ep < 100 keV (Band 
2003), and bursts with Ep < 100 keV arc easily lost. If it is 
really the case, the Ep distribution in our threshold setting 
might still suffer a little of biases for Ep < 100 keV. The 
upcoming Swift satellite, which will be scheduled for launch 
in 2004 September (Gehrels 2004), is marginally more sen- 
sitive than BATSE for Ep > 100 kcV but significantly more 
sensitive for Ep < 100 keV (Band 2003). It might establish 
a more reliable Ep distribution, which may be used for a 
further examination of our results. 

The true GRB rate as a function of rcdshift is difficult 
to determine from the current GRB sample. We assume that 
the GRB rate is proportional to the star formation rate. Var- 
ious models of star formation rates are presented in the liter- 
ature, and they arc quite different. In our analysis, we focus 
on the comparison of the simulated 9 distribution with that 
predicted by the US.] model. We find that the comparison 
result is not significantly affected by the redshift distribu- 
tion. In our analysis, the SFl model from Porciani & Madau 
(2001) is used. We use the SF2 model (Porciani & Madau 
2001) to perform the simulations. The results are shown in 
Figure 4. One can see that the comparison of the two distri- 
butions is not signiflcantly affected by the model of redshift 
distribution. 

Based on our results and the above discussion, we con- 
clude that both the probability distributions, P{0) and 
P{9,z), derived from our simulations are consistent with 
that predicted by the USJ model. This consistency may be 
regarded as a support of USJ model. 
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greatly. We also thank Bing Zhang and Xinyu Dai for their 



^ http: / /www.batse.msfc.nasa.gov/batse/ 




i9(rad) 



Figure 4. Same as Figure 1 but SF2 model is used. 
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